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Abstract 

Acute  stress  responsiveness  is  a  quantitative  trait  that  varies  in  severity  from  one  individual 
to  another;  however,  the  genetic  component  underlying  the  individual  variation  is  largely 
unknown.  Fischer  344  (F344)  and  Wistar  Kyoto  (WKY)  rat  strains  show  large  differences  in 
behavioral  responsiveness  to  acute  stress,  such  as  freezing  behavior  in  response  to  foot- 
shock  during  the  conditioning  phase  of  contextual  fear  conditioning  (CFC).  Quantitative  trait 
loci  (QTL)  have  been  identified  for  behavioral  responsiveness  to  acute  stress  in  the  defen¬ 
sive  burying  (DB)  and  open  field  test  (OFT)  from  a  reciprocal  F2  cross  of  F344  and  WKY  rat 
strains.  These  included  a  significant  QTL  on  chromosome  6  (StrespIO).  Here,  we  hypothe¬ 
sized  that  the  StrespIO  region  harbors  genes  with  sequence  variation(s)  that  contribute  to 
differences  in  multiple  behavioral  response  phenotypes  between  the  F344  and  WKY  rat 
strains.  To  test  this  hypothesis,  first  we  identified  differentially  expressed  genes  within  the 
StrespIO  QTL  in  the  hippocampus,  amygdala,  and  frontal  cortex  of  F344  and  WKY  male 
rats  using  genome-wide  microarray  analyses.  Genes  with  both  expression  differences  and 
non-synonymous  sequence  variations  in  their  coding  regions  were  considered  candidate 
quantitative  trait  genes  (QTGs).  As  a  proof-of-concept,  the  F344.WKY-Sfresp70congenic 
strain  was  generated  with  the  StrespIO  WKY  donor  region  into  the  F344  recipient  strain. 
This  congenic  strain  showed  behavioral  phenotypes  similar  to  those  of  WKYs.  Expression 
patterns  of  Gpatchll  (G-patch  domain  containing  11),  Cdkl4  (Cyclin  dependent  kinase  like 
4),  and  Drcl  (Dynein  regulatory  complex  subunit  1)  paralleled  that  of  WKY  in  the  F344. 
WKY-Sfresp70strain  matching  the  behavioral  profiles  of  WKY  as  opposed  to  F344  parental 
strains.  We  propose  that  these  genes  are  candidate  QTGs  for  behavioral  responsiveness  to 
acute  stress. 
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Introduction 

Acute  stress  is  a  common  occurrence  in  everybody’s  life  with  differing  severity.  Individual  var¬ 
iation  in  the  response  to  acute  stress  depends  on  the  sensitivity  to  the  stressor  and  the  ability  to 
cope  with  it.  Genetic  studies  support  the  assumption  that  acute  stress  responsiveness  is  a  quan¬ 
titative  trait  [1-3].  Its  narrow-sense  heritability  ranges  between  0.4-0.6  defined  by  multivariate 
genetic  analyses  in  family  and  twin  studies  [4-7].  The  heritability  of  acute  stress  responsiveness 
is  estimated  to  be  similar  in  rodents  to  that  in  humans  [8,  9].  However,  only  a  limited  number 
of  candidate  genes  have  been  proposed  for  individual  variation  in  behavioral  acute  stress 
responsiveness  to  date  [10-14].  The  purpose  of  this  study  was  to  identify  candidate  genes  with 
sequence  variations  that  contribute  to  variations  in  individual  stress  vulnerability. 

We  have  studied  the  genetic  components  of  behavioral  responsiveness  to  acute  stress  in  the 
reciprocally  crossed  F2  generation  of  Wistar  Kyoto  (WKY)  and  Fischer  344  (F344)  rat  strains 
using  the  defensive  burying  (DB)  and  open  field  tests  (OFT)  [15-19].  The  parental  WKY  strain 
consistently  exhibits  hypoactive  and  avoidant  behavior  compared  to  that  of  the  F344  [19].  We 
argue  that  this  consistent  repertoire  of  behavior  across  various  tests  is  not  composed  of  dis¬ 
crete,  individual  reactions  to  each  situation,  but  belies  a  more  global  response  to  an  acute  chal¬ 
lenge  with  substantial  genetic  origin.  Among  the  quantitative  trait  loci  (QTL)  found  in  these 
studies,  a  highly  significant  locus  on  chromosome  6,  StresplO,  has  been  identified  for  multiple 
phenotypes;  a  potential  pleiotropy  for  behavioral  stress  responsiveness  (Table  1)  [16].  The  phe¬ 
notypes  associated  with  this  QTL  include  latency  to  bury,  grooming,  rearing  in  DB,  and  rear¬ 
ing  in  OFT.  It  is  notable  that  a  discrepancy  was  described  for  the  OFT  and  DB  rearing  QTLs 
using  genetic  mapping  from  the  time  of  the  original  publication  [19].  However,  physical  map¬ 
ping  using  the  more  recent  rat  genome  version  Rnor_6.0  identified  that  these  QTLs  indeed 
overlap.  Furthermore,  rearing  measures  highly  correlate  between  the  two  tests  [19].  All  of 
these  behaviors  represent  a  general  stress  response  to  an  acute  stressor,  whether  it  be  the  shock 
in  the  DB  or  the  novelty  of  OFT.  Thus,  this  chromosomal  region  is  likely  to  harbor  one  or 
more  genes  with  sequence  variation(s)  that  contribute  to  phenotypic  variations  in  general 
stress  responsiveness  to  acute  stressors  between  these  strains. 

In  this  study,  we  aimed  to  identify  candidate  quantitative  trait  genes  (QTG)  that  contribute 
to  differences  in  acute  stress  responsiveness  between  F344  and  WKY  rat  strains  in  the  StresplO 
QTL.  To  aid  in  the  identification  of  QTG(s),  first  we  established  the  differentially  expressed 
genes  (DEGs)  between  the  two  parental  strains  within  the  StresplO  QTL  using  genome-wide 
microarray  analyses  in  the  hippocampus,  amygdala,  and  frontal  cortex.  Using  the  Rat  Genome 
Database  (www.rgd.mcw.edu)  and  the  original  sequences  of  the  two  parental  strains  [12, 14],  we 
then  determined  nonsynonymous  sequence  variations  in  the  coding  regions  of  DEGs  between 
the  F344  and  WKY  strains.  We  hypothesized  that  DEGs  that  show  expression  differences  paral¬ 
lel  to  behavioral  differences  between  the  two  strains  are  candidate  quantitative  trait  genes.  To 


Table  1.  List  of  behavioral  phenotypes  associated  with  the  StresplO  QTL. 


QTL  Symbol 

Position 

Behavior 

LOD 

Pointwise  significance 

Variance  (%) 

Reference 

StresplO 

6:1-35,623,029 

Latency  (DB) 

3.55 

0.00002 

3.5 

[15] 

Rearing  (DB)*A 

3.72 

0.00019 

1.6 

[19] 

Grooming  (DB)* 

3.1 

0.00079 

[19] 

Rearing  (OFT) 

8.24 

2.82E-10 

8.9 

[19] 

*  suggestive,  trait  x  lineage  covariates 
A  Rearing  (DB)  QTL  was  originally  identified  at  62cM 
From  the  Rat  Genome  Database  version  Rnor_6.0 

https://doi.org/10.1371/journal.pone.0194293.t001 
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investigate  the  generalizability  of  the  findings  to  a  different  acute  stress  phenotype,  we  measured 
the  behavioral  response  differences  in  the  conditioning  part  of  the  contextual  fear  conditioning 
(CFC)  paradigm  between  the  two  strains  since  the  acute  stress  of  the  electric  shock  in  the  CFC  is 
similar  to  that  shock  in  the  DB  test.  For  proof-of-concept,  we  measured  the  transcript  levels  of 
these  candidate  QTGs  in  the  brain  regions  of  a  congenic  strain  generated  by  inserting  the 
StresplO  WKY  donor  region  into  the  F344  recipient  strain.  DEGs  with  expression  profiles  in  the 
FM4.WKY -StresplO  congenic  strain  similar  to  the  WKY’s,  but  different  from  the  F344,  suggest 
that  the  sequence  variations  within  these  genes  might  contribute  to  the  StresplO  QTL. 

Materials  and  methods 

Animal  care  and  treatment 

All  animal  experiments  were  performed  according  to  procedures  approved  by  the  Institutional 
Animal  Care  and  Use  Committee  of  Northwestern  University.  Adult  male  Fischer  344  (F344) 
and  Wistar  Kyoto  (WKY)  rats  (approximately  3  months  old)  were  obtained  from  Envigo 
(Indianapolis,  IN)  and  Charles  River  Laboratories  (Wilmington,  MA),  respectively.  All  animals 
were  group-housed  (2-3  per  cage)  in  a  temperature-controlled  environment  with  12  h  light/ 
dark  cycles  and  allowed  feed  and  water  ad-libitum. 

The  experimental  design  required  three  animal  cohorts.  Specifically,  one  cohort  of  control 
F344  and  WKY  rats  were  not  subjected  to  any  behavioral  testing.  RNA  was  isolated  from  the 
three  brain  regions  and  used  for  the  microarray  experiments.  The  second  cohort  consisted  of 
F344  and  WKY  rats  that  were  subjected  to  CFC  testing  and  their  RNA  was  isolated  from  the 
three  brain  regions.  This  RNA  was  used  for  expression  analyses  by  RT-PCR  of  candidate  genes 
together  with  RNA  isolated  from  the  1  • 3 44.  W  K  Y  -  .S’  tresp  1 0  rats.  The  third  cohort  of  F344, 

WKY,  and  F344.WKY-SfresplO  rats  were  maintained  independently  and  used  for  only  behav¬ 
ioral  testing  in  the  OFT  followed  by  DB  three  weeks  later. 

Construction  of  congenic  F344.WKY-StresplO 

The  F344.  WKY -Stresp  1 0  congenic  strain  was  generated  by  repeated  marker-directed  back- 
crossing  of  the  F344-WKY  F,  generation  into  the  F344  parental  strain.  The  markers  used  are 
described  before  [16, 19].  The  N15  generation  of  F344.WKY -Stresp  10  congenic  male  animals 
(approximately  3  months  old)  were  used  for  behavioral  experimentation  together  with  the  simul¬ 
taneously  maintained  F344  and  WKY  rats.  The  15th  generation  of  the  1  •’ 3 44 .  W K  Y  -  .S’ tresp  1 0  strain 
was  homozygous  from  135,053  to  6,709,713  bp  and  from  19,464,437  to  the  end  of  the  QTL  at 
28,931,796  bp,  and  heterozygous  from  6,963,239  to  15,597,330  bp  (SI  Table).  Additionally,  after 
15  generations  of  backcross,  congenic  strains  are  known  to  have  less  then  5%  heterozygosity 
throughout  the  genome  [20].  This  strain  was  characterized  using  the  10K  Affymetrix  Targeted 
Genotyping  Array  (Affymetrix,  Santa  Clara,  CA)  [21].  Using  this  array,  4  pg  of  rat  DNA  was  gen- 
otyped  with  the  GeneChip  Scanner  2000  Targeted  Genotyping  System  (GTGS)  following  the 
manufacturer’s  protocol  (Affymetrix,  Santa  Clara,  CA)  [22].  The  genotyping  data  was  analyzed 
using  the  corresponding  software  for  GTGS  (Affymetrix,  Santa  Clara,  CA).  The  10K  Affymetrix 
Targeted  Genotyping  Array  used  the  rat  genome  version  RGSC_3.4.  The  QTL  locations  were 
translated  to  the  latest  version  or  the  rat  genome,  Rnor_6.0.  SI  Table  lists  the  QTL  locations 
from  both  versions  of  the  rat  genome. 

Behavioral  tests 

Contextual  fear  conditioning.  The  fear  conditioning  phase  of  the  contextual  fear  condi¬ 
tioning  test  was  carried  out  as  described  previously  [23].  Briefly,  male  rats  of  the  F344  and 
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WKY  strains  were  individually  placed  in  the  CFC  chamber  (Technical  &  Scientific  Equipment, 
Bad  Homburg,  Germany)  and  exposed  to  3  minutes  of  habituation  followed  by  3  footshocks 
(0.8  mA,  1  sec  duration,  spaced  1  minute  apart).  The  behavior  observed  after  the  footshocks 
was  analyzed  by  the  TSE  Videomot  software  version  5.75  (Bad  Homburg,  Germany),  which 
measures  the  animals’  movement,  including  distance  travelled,  number  of  times  rearing,  and 
freezing  duration.  The  rats  were  sacrificed  by  quick  decapitation  24  hours  later.  Whole  brains 
were  collected  in  RNAlater  (Ambion,  Austin,  TX)  and  frozen  at  -80°C  until  dissection  of  brain 
regions.  Total  RNA  was  isolated  from  the  hippocampus,  amygdala,  and  frontal  cortex,  and  was 
used  for  the  expression  analyses  of  candidate  genes.  All  rats  were  approximately  3  months  old 
at  the  time  of  sacrifice  and  tissue  sample  collection. 

Open  field  test.  The  open  field  test  was  carried  out  on  another  cohort  of  animals,  as  previ¬ 
ously  described  [23].  Briefly,  male  rats  of  the  F344,  WKY,  and  F344.WKY -StresplO  strains 
were  placed  in  a  circular  arena  (82  cm  in  diameter)  with  a  30  cm  high  wall,  lit  to  a  brightness 
of  approximately  60  lux  by  indirect  overhead  lighting.  The  arena  contained  an  inner  concen¬ 
tric  circle  (50  cm  in  diameter)  designated  as  the  inner  zone.  The  rats  were  placed  in  the  center 
of  the  arena  and  allowed  to  move  freely  for  10  min  while  the  activity  was  recorded  and  tracked 
by  the  TSE  Videomot  software  2  version  5.75  (Bad  Homburg,  Germany),  which  measured  the 
number  of  times  rearing  by  the  animal. 

Defensive  burying.  Three  weeks  after  the  OFT,  the  defensive  burying  test  was  carried  out 
on  the  same  animals  as  previously  described  in  the  QTL  studies  [15,  16, 19].  Briefly,  male  rats 
of  the  three  strains  were  habituated  to  a  Plexiglas  chamber  (40  cm  square,  60  cm  high)  with 
bedding  (wood  shaving)  (7  cm  deep,  1  cm  below  the  hole  for  the  prod)  for  15  min  each  day, 
for  three  consecutive  days,  between  10:00  AM  and  2:00  PM.  On  the  fourth  day,  a  continuously 
electrified  prod  was  introduced  into  the  chamber,  which  delivered  a  4.5  mA  shock  when  the 
rat  touched  it.  The  shock  was  generated  from  a  shock  generator  (Lafayette  Instruments,  San 
Diego,  CA)  set  at  4.5  mA.  The  rats  typically  explored  the  novel  prod  and  received  a  shock, 
which  started  the  15-minute  video-taped  test.  Once  shocked,  animals  typically  retreated  to  the 
back  of  the  cage,  and  either  remained  there  (classic  WKY  behavior)  or  began  spraying  bedding 
toward  the  prod  in  an  effort  to  cover  it.  Behaviors  scored  by  an  observer  blind  to  the  genotype 
of  the  animal  include  the  latency  to  begin  burying,  the  total  time  spent  burying  (duration  of 
burying),  the  duration  of  grooming,  and  the  number  of  times  rearing  by  the  animal.  The  rats 
were  sacrificed  by  quick  decapitation  two  weeks  after  the  test.  Whole  brains  were  collected  in 
RNAlater  (Ambion,  Austin,  TX)  and  frozen  at  -80”C  until  dissection.  All  rats  were  approxi¬ 
mately  3  months  old  at  the  time  of  sacrifice  and  tissue  sample  collection. 

Brain  dissection  and  RNA  isolation 

Brain  regions  (hippocampus,  amygdala,  and  frontal  cortex)  were  dissected  from  adult  male 
F344,  WKY  and  F344.  WKY-Sfresp  1 0  congenic  rats  as  described  previously  [24]  and  stored  in 
RNAlater  (Ambion,  Austin,  TX)  at  -80  °C.  Briefly,  a  rat  brain  matrix  was  used  to  prepare  pri¬ 
mary  sections  from  which  the  brain  regions  were  dissected  according  the  Paxinos  coordinates 
[25]:  hippocampus  (AP-2.12  to  -6.0,  ML  0  to  5.0,  DV  5.4  to  7.6),  amygdala  (AP  -0.58  to  -2.18, 
ML  1.5  to  4.5,  DV  4  to  5.75),  and  frontal  cortex  (AP  5.20  to  1.70,  ML  0  to  3.3,  DV  9.0  to  4.4). 

Total  RNA  was  isolated  from  individual  tissue  samples  using  the  Direct-zol  RNA  MiniPrep 
Kit  (Zymo  Research,  Irvine,  CA).  Briefly,  individual  tissue  samples  were  homogenized  in  TRI- 
zol  reagent  (Ambion,  Austin,  TX)  and  RNA  was  immediately  isolated  using  the  kit  following 
the  manufacturer’s  protocol.  Total  RNA  (2  pg)  was  reverse  transcribed  to  generate  double- 
stranded  cDNA  using  the  Superscript  VILO  cDNA  Synthesis  Kit  (ThermoFisher  Scientific, 
Waltham,  MA). 
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Microarray  analysis 

Genome-wide  microarray  analysis  was  performed  from  RNA  isolated  from  hippocampi, 
amygdalae,  and  frontal  cortexes  of  F344  and  WKY  adult  male  rats  that  were  unstressed,  with¬ 
out  undergoing  any  behavioral  testing,  as  described  previously  [26].  Briefly,  the  cDNA  gener¬ 
ated  from  the  RNA  was  linearly  amplified  and  labeled  with  biotinylated  nucleotides  in  an  in 
vitro  transcription  reaction  using  the  Illumina  TotalPrep  RNA  Amplification  Kit  (San  Diego, 
CA)  to  make  cRNA.  1.5  pg  of  biotin-labelled  and  fragmented  cRNA  was  then  hybridized  onto 
RatRef-12  Expression  BeadChips  (Illumina,  San  Diego,  CA).  The  BeadChips  have  multiple 
probes  per  transcript  to  mitigate  hybridization  bias  artifacts.  Probe  intensity  data  from  the 
BeadChips  were  directly  read  into  the  R  software  environment  (http://www.R-project.org) 
from  bead  summary  files  produced  by  BeadStudio  using  the  R/beadarray  package  [26,  27]. 
Quantile  normalization  was  applied  to  the  data  using  the  R/preprocessCore  package  [26,  28]. 
Data  quality  was  assessed  using  histograms  of  signal  intensities,  scatter  plots,  and  hierarchical 
clustering  of  samples,  as  previously  described  [26].  Statistical  significance  of  microarray 
expression  differences  between  F344  and  WKY  was  determined  using  ANOVA  methods 
within  the  R/maanova  package  as  previously  described  [26,  29].  DEGs  were  determined 
between  strains  with  an  FDR-adjusted  P-value  less  than  0.05  and  a  fold  change  greater  than 
1.3  (30%  increase  or  decrease).  This  criterion  has  been  well-established  to  give  biologically 
meaningful  datasets  when  interpreting  differential  gene  expression  profiles  in  microarray 
experiments  [30,  31].  The  microarray  data  used  the  rat  genome  RGSC_v3.4  for  identification 
of  transcripts,  which  was  translated  to  the  latest  version  of  rat  genome,  Rnor_6.0. 

The  differential  gene  expression  profiles  between  the  two  strains  were  determined  with  a 
significance  criteria  of  an  FDR-adjusted  P  value  less  than  0.05  and  an  absolute  fold  change 
above  1.3  (30%  increase  or  decrease).  With  these  criteria,  we  found  1,030  DEGs  in  the  hippo¬ 
campus,  769  in  the  amygdala,  and  976  in  the  frontal  cortex,  as  listed  in  S3-S5  Tables. 

Identification  of  single  nucleotide  polymorphisms  in  coding  regions 

The  genes  with  non-synonymous  single  nucleotide  polymorphisms  (SNPs)  within  coding 
regions  between  the  F344  and  WKY  genomes  were  obtained  from  the  Rat  Genome  Sequencing 
and  Mapping  Consortium  and  Baud  et  al.  [12, 14].  In  these  studies,  both  F344  and  WKY 
genomes  were  first  mapped  to  the  Brown  Norway  reference  genome,  version  RGSC_3.4. 

Using  the  Integrative  Genomics  Viewer  (Broad  Institute,  Cambridge,  Massachusetts),  we  iden¬ 
tified  the  coding  sequence  variations  between  the  F344  and  WKY  genomes  and  translated 
them  to  the  rat  genome  version  Rnor_6.0.  For  quality  control,  we  set  the  criteria  that  the  cover¬ 
age  for  each  SNPs  had  to  be  greater  than  10  reads,  where  reads  are  used  to  reconstruct  the 
sequence.  The  more  reads  a  sequence  had,  the  more  reliable  the  data.  Furthermore,  the  single 
nucleotide  variations  had  to  be  called  in  at  least  50%  of  the  reads. 

Among  inherited  gene  variations  in  humans,  nonsynonymous  single  nucleotide  polymor¬ 
phisms  that  lead  to  an  amino  acid  change  in  the  protein  product  are  most  relevant  to  human 
inherited  diseases  [32].  Therefore,  as  the  first  step  in  the  identification  of  candidate  QTGs,  we 
focused  on  this  subset  of  genes.  Genes  with  non-synonymous  SNPs  within  the  heterozygous 
StresplO  region  were  excluded  from  our  analysis.  Of  the  16  DEGs  mapped  within  the  StresplO 
region,  10  genes  were  found  to  contain  SNPs  (Table  2).  SNPs  are  listed  in  S9  Table. 

Real-time  RT-PCR 

The  real-time  reverse  transcription-polymerase  chain  reaction  (qRT-PCR)  was  carried  out 
using  RNA  from  the  second  cohort  of  F344  and  WKY  rats  that  were  subjected  to  CFC  testing, 
as  well  as  RNA  from  the  F344.WKY-Sfresp  1 0  rats.  Primers  for  each  gene  were  designed  using 
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Table  2.  Differentially  expressed  genes  within  the  StresplO  QTL  region  containing  coding  region  non-synonymous  SNPs  between  F344  and  WKY. 


Gene  Symbol 

Gene  Description 

Chr 

Start 

End 

P- value 

0  =  <  1.00E-08 

Gpatchll^’^ 

G-patch  domain  containing  1 1 

6 

1410507 

1423041 

h:  0,  a:  0 ,/:  0 

PrkdV 

Protein  kinase  D3 

6 

1546018 

1622232 

f:  1.18E-07 

Cyplbl^ 

Cytochrome  P450,  family  1,  subfamily  b,  polypeptide  1 

6 

2307808 

2316722 

h:  0,/:  0 

Cdkl4'’ 

Cyclin-dependent  kinase-like  4 

6 

3234090 

3254779 

h:  7.39E-05 

Mta3Kf 

Metastasis  associated  1  family,  member  3 

6 

6908684 

7031828 

h  :  0 ,/:  0 

Ttc27“'f 

Tetratricopeptide  repeat  domain  27 

6 

21735834 

21880003 

a:  1.42E-07,/:  0 

A5khAf 

Anaplastic  lymphoma  receptor 

6 

22696397 

23203775 

h :  0,  a:  0,/:  0 

Rbks'1 

Ribokinase 

6 

26051396 

26128906 

h:  0 

Drcl0 

Dynein  regulatory  complex  subunit  1 

6 

27425235 

27460038 

a:  4.74E-08 

Dtnb,,'“/ 

Dystrobrevin,  beta 

6 

27975417 

28177214 

h:  5.09E-07,  a:  3.79E-05,/i  1.56E-05 

h:  hippocampus,  a:  amygdala,  /’:  frontal  cortex 
From  the  Rat  Genome  Database  version  Rnor_6.0 

https://doi.org/10.1371/journal.pone.0194293.t002 


Primer  Express  Software  version  3.0  (Applied  Biosystems,  Carlsbad,  CA).  The  default  setting 
was  used  to  design  primers  that  amplify  80  to  150  bp  regions.  The  primer  sequences  for  the 
candidate  quantitative  trait  genes  are  listed  in  S2  Table.  Five  ng  of  cDNA  were  amplified  in 
20  pL  reactions  (IX  SYBR  Green  Master  Mix  (ThermoFisher  Scientific),  250  pM  primers)  in 
QuantStudio  7  Flex  Real-Time  PCR  System  (ThermoFisher  Scientific,  Waltham,  MA)  using 
the  relative  quantification  (-AACt)  method,  with  Gapdh  as  the  housekeeping  gene  and  a  gen¬ 
eral  calibrator.  We  have  established  that  there  was  no  change  in  hippocampal  Gapdh  expres¬ 
sion  across  strains  and  conditions. 

Statistical  analysis 

All  data  were  presented  as  mean  ±  standard  error  of  mean.  Outliers  from  the  quantitative  RT- 
PCR  data  were  determined  as  being  more  than  two  standard  deviation  away  from  the  mean. 
Therefore,  the  number  of  samples  per  group  differ  between  target  genes.  All  statistical  analyses 
were  performed  using  GraphPad  Prism  v  7.0  (GraphPad  Software,  La  Jolla,  CA).  Statistical  signifi¬ 
cance  of  differences  between  strains  were  determined  by  ANOVA,  followed  by  post-hoc  analysis 
with  the  Bonferroni’s  correction  for  multiple  comparisons.  Statistical  significance  was  considered 
at  an  adjusted  P-  value  of  less  than  0.05.  When  significant  main  effects  were  indicated  by  the 
ANOVA,  but  the  Bonferroni’s  multiple  comparisons  test  did  not  show  significance,  hypothesis 
testing  by  Student’s  /-test  was  carried  out  between  groups.  Our  decision  to  apply  the  Student’s  t- 
test  was  based  on  an  increasing  number  of  discussions  arguing  that  P-values  are  not  as  reliable  as 
it  is  previously  thought  [33]  and  that  while  a  three-group  comparison  ANOVA  may  not  result  in 
significance,  two  groups  of  the  three  can  differ  from  each  other  at  the  P  <  0.05  level  [34], 

ANOVA  results  are  given  in  the  results  and  post-hoc  significances  are  noted  in  the  figures. 

Results 

Microarray  analysis 

Genome-wide  microarrays  were  performed  on  hippocampal,  amygdalar,  and  frontal  cortex 
RNA  from  unstressed  F344  and  WKY  male  rats.  To  identify  potential  QTGs,  DEGs  were 
mapped  within  the  StresplO  QTL  chromosomal  region.  This  region  is  associated  with  multiple 
behavioral  phenotypes  in  response  to  acute  stress  (Table  1).  The  chromosomal  location  of 
these  QTGs  was  mapped  from  RGSC_v3.4  to  the  latest  version  of  rat  genome  Rnor_6.0.  Of  the 


PLOS  ONE  |  https://doi.org/10.1371/journal.pone.0194293  March  12,  2018  6/18 

Distribution  A.  Approved  for  public  release  (PA):  distribution  unlimited. 


ONE 


Candidate  quantitative  trait  genes  for  acute  stress  responsiveness 


>PLOS 


DEGs,  14  genes  in  the  hippocampus,  12  in  the  amygdala,  and  14  in  the  frontal  cortex  were 
found  within  the  StresplO  region  (S6-S8  Tables).  Since  many  of  these  genes  overlap  in  two  or 
more  brain  regions,  there  were  a  total  of  18  genes  to  investigate.  Precisely,  there  were  six  genes 
( Gpatchll ,  Slc3al,  Camkmt,  Aik,  Dtnb,  Klhl29 )  that  overlapped  in  all  three  brain  regions: 
Rasgrp3  between  the  hippocampus  and  amygdala,  Cypl  bl  and  Mta3  between  the  hippocampus 
and  frontal  cortex,  and  Ttc27  between  the  amygdala  and  frontal  cortex. 

Single  nucleotide  polymorphisms  in  coding  regions  of  DEGs  mapped 
within  Strespl  0 

To  determine  which  genes  contribute  to  the  phenotypic  variation  in  acute  stress  responsiveness, 
we  identified  the  genes  with  non-synonymous  single  nucleotide  polymorphisms  (SNPs)  within 
coding  regions  between  the  F344  and  WKY  genomes,  which  were  obtained  from  Baud  et  al.  [12, 
14],  The  SNPs  were  identified  first  between  the  WKY,  F344  and  the  Brown  Norway  reference 
genome  and  then  between  the  WKY  and  F344  genome  using  Rnor_6.0.  The  candidate  SNPs  had 
to  be  greater  than  10  number  of  reads  and  the  single  nucleotide  variations  in  greater  than  50%  of 
the  reads.  Of  the  16  DEGs  mapped  within  the  homozygous  StresplO  region,  10  genes  were  found 
to  contain  non-synonymous  SNPs  in  coding  regions  (Table  2).  These  SNPs  are  listed  in  S9  Table. 

Acute  stress  responsiveness  in  F344  and  WKY  rat  strains 

To  investigate  the  generalizability  of  the  findings  to  a  different  acute  stress  phenotype,  we  mea¬ 
sured  the  acute  stress  response  differences  between  the  two  strains  in  the  conditioning  part 
(day  1)  of  the  CFC  (Fig  1).  Similar  to  behaviors  in  the  DB  and  OFT  [15,  16, 19],  WKY  rats 
exhibited  significantly  more  freezing  in  response  to  the  stress  of  the  footshock  (f(30)  =  4.066; 

P  <  0.01),  which  is  a  clear  indication  of  a  more  passive  stress  response  (Fig  1A).  Furthermore, 
WKYs  also  exhibited  a  significantly  lower  frequency  of  rearing  (1(29)  =  4.00;  P  <  0.01),  which 
is  an  avoidance  response  [19,  35].  WKYs  also  exhibited  a  more  hypoactive  response  to  the  foot- 
shocks,  measured  by  significantly  shorter  distance  travelled  in  the  chamber  compared  to  the 
F344  rats  (f(30)  =  1.813;  P  <  0.01)  (Fig  IB  and  1C). 


Confirmation  of  Strespl  0  phenotypes  in  F344  and  WKY  parental  strains, 
and  in  the  F344.WKY-Sfre.spf  Ocongenic  strain 

The  phenotypes  associated  with  StesplO  include  latency  to  bury,  grooming,  and  rearing  in  the 
DB  test,  and  rearing  in  the  OFT;  all  of  which  represent  an  acute  stress  response.  To  confirm 
that  the  StresplO  QTL  contributes  to  the  variation  in  these  behavioral  acute  stress  responses,  a 
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Fig  1.  Behavioral  responses  to  footshock  of  F344  and  WKY  adult  male  rats.  (A)  WKY  adult  male  rats  froze  longer  (freezing  duration  in  seconds)  following  the 
footshocks  in  the  CFC  conditioning  chamber.  (B)  Number  of  rears  are  significantly  lower  for  WKYs  compared  to  F344s.  (C)  F344s  move  around  the  chamber  after 
the  footshock  significantly  more  than  WKYs  measured  by  distance  traveled  in  cm.  Data  are  presented  as  mean  ±  standard  error  of  mean.  N  =  16  per  group.  ** 

P  <  0.01  by  Student’s  f-test. 


https://doi.org/10.1371/journal.pone.0194293.g001 
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congenic  strain,  F344.WKY-SfresplO,  was  generated  (SI  Table).  Fig  2  shows  the  phenotypic 
differences  between  the  F344,  WKY,  and  F344.WKY-Sfrespl  0  strains.  In  the  DB  test,  both  the 
WKY  and  F344.WKY-Sfresp  1 0  rats  exhibited  enhanced  avoidance  and  hypoactive  responses  to 
acute  stress,  measured  by  significantly  longer  latency  to  bury  (F[3,46]  =  44.55;  P  <  0.01)  and 
shorter  duration  of  burying  the  electrified  prod  fF[2,43]  =  23.96;  P  <  0.01)  compared  to  the 
F344s  (Fig  2A  and  2B).  Furthermore,  WKY  and  F344. WKY -Stresp  1 0  rats  were  less  active 
showing  less  grooming  (strain:  F[2,43]  =  11.80,  P  <  0.01)  and  rearing  (strain:  .F[2,25]  =  34.04, 

P  <  0.01)  compared  to  the  F344  rats  (Fig  2C  and  2D).  In  the  OFT,  WKYs  reared  significantly 
less  (F[2,87]  =  14.79,  P  <  0.01)  and  the  number  of  rearing  for  the  F344.WKY-SfresplO  were 
intermediary  not  differing  from  either  parental  strain  significantly  (Fig  2E). 


Quantitative  expression  of  DEGs  with  sequence  variations  in  the  brain 
regions  of  F344,  WKY,  and  F344-WKY. Stresp  1 0  adult  male  rats 

Quantitative  expression  analyses  of  the  10  genes  (Table  2)  were  carried  out  using  RNA  from 
the  hippocampus,  amygdala,  and  frontal  cortex  of  F344,  WKY,  and  congenic  strain  F344. 
WKY -Stresp  10  male  rats. 


Latency  to  Bury  in  DB  Duration  of  Burying  in  DB  (^)  Grooming  in  DB 


Fig  2.  Phenotypic  differences  between  the  F344,  WKY,  and  F344.WKY-Sfrespl0  strains  in  the  defensive  burying  and  open  field  tests.  (A,  B)  Both  WKY  and  F344. 
WKY -Stresp  10  adult  male  rats  exhibited  significantly  longer  latency  to  bury  and  shorter  duration  of  burying  the  electrified  prod  compared  to  the  F344s.  (C,  D)  Both 
WKY  and  F344.WKY -Stresp  10  rats  groomed  less  (time  spent  grooming  in  seconds)  and  reared  less  compared  to  the  F344s  in  the  DB  paradigm.  (E)  WKYs  also  rear 
significantly  less  than  F344s  in  the  OFT,  while  the  number  of  rearing  of  F344.WKY-Sfrespl0  did  not  differ  significantly  from  either  parental  strain.  Data  are  presented  as 
mean  ±  standard  error  of  mean.  N  =  13-18  per  group.  **  P  <  0.01  by  two-way  ANOVA  followed  by  Bonferroni’s  multiple  comparisons  test. 

https://doi.org/10.1371/journal.pone.0194293.g002 


PLOS  ONE  |  https://doi.org/10.1371/journal.pone.0194293  March  12,  2018  8/18 

Distribution  A.  Approved  for  public  release  (PA):  distribution  unlimited. 


ONE 


Candidate  quantitative  trait  genes  for  acute  stress  responsiveness 


PLOS 


(B)  Cyplbl 


Fig  3.  Expression  analyses  of  candidate  quantitative  trait  genes  in  the  hippocampus,  amygdala,  and  frontal  cortex  of  F344,  WKY,  and  F344.WKY-StresplO  adult 
male  rats.  (A)  Candidate  genes  with  similar  expression  between  the  WKY  and  F344.WKY-StresplO  strains,  but  different  from  the  F344,  in  two  or  more  brain  regions. 
(B)  Genes  that  express  a  similar  pattern  in  only  one  brain  region.  Data  are  presented  as  mean  ±  standard  error  of  mean.  N  =  5-7  per  group.  *  P  <  0.05  and  **  P  <  0.01 
by  two-way  ANOVA  followed  by  Bonferroni’s  multiple  comparisons  test.  #  P  <  0.05  and  ##  P  <  0.01  by  Student’s  f-test. 


https://doi.org/10.1371/journal.pone.0194293.g003 


Genes  with  similar  expression  between  the  WKY  and  F344.WKY -StresplO  strains,  but  dif¬ 
ferent  from  the  F344,  in  two  or  more  brain  regions  were  considered  to  be  strong  candidate 
QTGs  that  likely  contribute  to  the  variation  in  acute  stress  responsiveness  between  the  parental 
strains.  These  candidate  QTGs  include  Gpatchll  (G-patch  domain  containing  11),  Cdkl4 
(Cyclin  dependent  kinase  like  1),  and  Drcl  (Dynein  regulatory  complex  subunit  1)  (Fig  3A). 
Expression  of  Gpatchll  differed  significantly  by  brain  region  and  strain  (brain  region:  F[2,51] 

=  12.75,  P  <  0.01;  strain:  F[2,51]  =  21.82,  P  <  0.01).  Expression  of  Drc  differed  significantly 
by  brain  region  and  strain  as  well  the  interaction  between  the  two  (brain  region:  F[2,43]  = 
15.39,  P  <  0.01;  strain:  F[2,43]  =  96.56,  P  <  0.01;  interaction:  F[4,43]  =  3.844,  P  <  0.01).  Both 
Gpatchll  and  Drcl  showed  significant  differences  in  expression  between  the  F344s  and  WKYs 
and  F344s  and  1-344.  WKY -Stresp  1 0  rats  in  all  three  brain  regions.  Expression  of  Cdlk4  differed 
significantly  by  brain  region  and  strain  as  well  as  the  interaction  between  the  two  (brain  re¬ 
gion:  F[2,44]  =  22.85,  P  <  0.01;  strain:  E[2,44]  =  21.82,  P  <  0.01,  interaction:  F[4,44]  =  3.34, 

P  <  0.05).  Cdkl4  also  showed  the  same  pattern  of  expression  in  the  hippocampus  and  amyg¬ 
dala,  while  transcript  levels  in  the  frontal  cortex  were  relatively  low. 

Additionally,  genes  that  express  the  same  pattern  but  only  in  one  brain  region  were  consid¬ 
ered.  Sequence  variations  in  these  genes  can  interfere  with  gene  expression  by  interacting  with 
brain  region-specific  modulators,  such  as  miRNAs,  transcription  factors,  and  binding  proteins, 
present  or  absent  in  the  specific  brain  region,  which  may  or  may  not  induce  transcription  [36, 
37].  Only  two  genes,  Cyplbl  (Cytochrome  P450  family  1  subfamily  B  member  1)  and  Rbks 
(ribokinase),  showed  this  pattern  of  expression.  (Fig  3B).  Brain  region-  and  strain-specific 
expression  of  Cyplbl  showed  a  complex  pattern  (brain  region:  F[2,45]  =  26.01,  P  <  0.01;  strain: 
F[2,45]  =  7.465,  P  <  0.01;  interaction:  E[4,45]  =  3.567,  P  <  0.05).  Cyplbl  expression  was  signifi¬ 
cantly  greater  in  the  WKY  and  the  F344.WKY-Sfresp  1 0  strains  compared  to  F344,  but  only  in 
the  hippocampus.  Similar  to  the  expression  of  Cyplbl,  Rbks  transcript  levels  were  regulated 
in  a  brain  region-  and  strain-specific  manner  (brain  region:  F[2,50]  =  48.67,  P  <  0.01;  strain: 
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-F[2,50]  =  10.43,  P  <  0.01;  interaction:  _F[4,50]  =  4.041,  P  <  0.01).  Rbks  was  expressed  at  signifi¬ 
cantly  lower  levels  in  the  WKY  and  1;344. VVKY -Stresp  1 0  rats  relative  to  the  F344  rats  in  the  fron¬ 
tal  cortex;  the  hippocampal  expression  profile  was  similar  but  F344  and  the  F344.WKY-SfresplO 
showed  significance  only  as  hypothesis  testing  (f(ll)  =  2.57,  P  <  0.01).  There  were  no  signifi¬ 
cant  differences  in  the  expression  of  this  gene  in  the  amygdala.  Prkd3,  Mta3,  Ttc27,  Aik,  and 
Dtnb  showed  no  strain-specific  effects. 


Discussion 

The  major  findings  of  this  study  point  to  candidate  QTGs  that  can  contribute  to  differences  in 
behavioral  responsiveness  to  acute  stress  between  the  F344  and  WKY  strains  across  multiple 
paradigms  and  phenotypes.  Brain  region-specific  gene  expression  differences  between  the  F344 
and  WKY  strains  and  a  congenic  strain,  incorporating  a  WKY  donor  chromosomal  regions 
mapped  to  multiple  stress  responsive  QTLs  into  the  F344  recipient  background,  identified  can¬ 
didate  genes  within  this  locus.  Those  genes  that  showed  sequence  variations  between  the  two 
parental  strains  are  proposed  to  be  candidate  QTGs.  These  candidate  QTGs  were  identified 
using  a  multifaceted  approach  that,  to  our  knowledge,  has  not  been  used  previously.  Specifically, 
this  approach  included  genome-wide  microarray  analyses  to  identify  DEGs  within  the  QTL  in 
question,  followed  by  non-synonymous  SNPs  within  these  DEGs,  the  generation  of  the  con¬ 
genic  strain  for  the  behavioral  consequences  associated  with  the  specific  QTL,  and  finally  the 
expression  analysis  involving  all  three  strains  with  a  different  analytical  method,  the  qPCR. 

Human  studies  have  identified  few  genetic  variations  associated  with  individual  differences 
in  behavioral  responses  to  acute  stress  (for  example,  [38,  39]).  All  of  these  candidate  genetic 
variations  were  associated  with  responses  to  fearful  faces  measured  by  multiple  imaging  and 
other  methods.  In  contrast,  animal  studies  have  identified  multiple  QTL  for  behavioral  respon¬ 
siveness  to  acute  stress,  but  very  few  QTGs  have  been  proposed  [11-14,  40,  41].  The  usual 
methods  to  detect  QTL,  such  as  backcrosses,  F2  crosses,  and  consomic  strains  [42-44],  usually 
identify  large  genomic  regions,  with  large  number  of  genes  mapped  to  them.  Other  techniques 
such  as  recombinant  inbred  strains,  congenic  strains,  and  heterogeneous  stocks  [10,  45]  are 
able  to  identify  smaller  genomic  regions,  but  still  have  no  confirmed  QTGs  for  acute  stress 
responsiveness. 

In  our  previous  QTL  studies,  using  the  recombinant  F2  generation  of  the  reciprocally 
crossed  F344  and  WKY,  we  identified  multiple  QTL  for  behavioral  responsiveness  to  acute 
stress  in  different  paradigms.  Specifically,  the  chromosomal  region  of  the  Stresp  10  QTL  was 
associated  with  latency  to  bury,  grooming,  and  rearing  in  the  DB  test,  and  rearing  in  the  OFT 
from  this  cross.  All  of  these  phenotypes  describe  a  behavioral  response  of  the  animal  either  to  a 
novel  environment  or  to  an  aversive  shock  stimulus,  which  are  characteristically  either  active 
or  passive.  In  the  DB  test,  these  options  lead  to  the  active  behaviors  of  avoiding  the  shock  by 
increasing  the  duration  of  burying  and  of  rearing,  or  the  passive  response  to  the  shock  mani¬ 
fested  by  freezing  and  thereby  increased  latency  to  bury  [46,  47],  The  acute  stress  in  the  OFT  is 
the  novel  environment  from  which  the  animal  is  trying  to  escape  by  rearing;  an  active  response 
to  this  stress.  These  overlapping  QTL  within  the  StresplO  locus  may  represent  pleiotropy,  or 
multiple  sequence  variations  interacting,  resulting  in  a  common  genetic  architecture  underly¬ 
ing  different  behavioral  responses  to  acute  stress. 

In  the  fear  conditioning  component  of  the  CFC  test,  the  animal  can  either  freeze  (parallel  to 
latency  to  bury  in  the  DB),  explore  (parallel  to  burying),  or  rear,  after  it  receives  the  footshock. 
Exploration  is  an  active  behavior  that  is  inversely  related  to  freezing,  while  rearing  is  a  risk- 
assessing  behavior  that  encompasses  exploratory,  activity,  and  excitability  components  [15, 16, 
19].  In  both  the  CFC  and  the  DB  tests,  the  initial  period  of  exploration  allows  the  animal  to 
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form  an  associative  memory  between  the  context  and  the  conditioning  stimulus.  This  associa¬ 
tion  has  been  proven  to  be  necessary  to  trigger  conditioned  fear  responses  in  the  CFC  and  is 
essential  to  the  learning  component  of  the  fear  conditioning  paradigm  [48,  49].  Thus,  the  indi¬ 
vidual,  genetic  differences  in  these  responses  likely  influence  the  degree  of  fear  memory  recall 
and,  thereby,  are  of  major  significance.  Regarding  the  DB  test,  re-exposure  to  the  DB  chamber 
without  the  shock,  similar  to  CFC,  suggested  that  this  is  also  a  learning  paradigm  [50].  Al¬ 
though  the  CFC  phenotypes  may  not  be  mapped  to  StresplO,  considering  the  abovementioned 
parallels  and  the  third  rearing  phenotype  showing  similar  differences  between  the  strains,  we 
presumed  that  these  acute  stress  phenotypes  are  relevant  to  the  current  study.  Indeed,  the  pat¬ 
terns  of  candidate  gene  expressions  reflect  the  behavioral  measures  in  the  CFC,  OFT,  and  DB 
tests,  suggesting  that  the  candidate  QTGs  contribute  to  variations  in  the  general  behavioral 
response  to  acute  stress,  as  we  had  hypothesized. 

In  this  study,  we  aimed  to  identify  candidate  QTGs  within  the  StresplO  QTL.  We  mined 
our  previously  collected  genome-wide  transcriptomic  data  to  find  genes  in  this  region  with 
brain  region-specific  expression  that  parallel  the  behavioral  responses  to  acute  stress  in  the 
F344  and  WKY  strains.  The  brain  regions  explored  are  the  hippocampus,  amygdala,  and  fron¬ 
tal  cortex,  which  are  all  intimately  involved  in  the  behaviors  discussed.  The  neural  circuitry 
connecting  the  amygdala  and  frontal  cortex  are  involved  in  the  emotional  responses  to  acute 
stress  [51-54].  The  ventral  hippocampus  is  involved  in  anxiety-like  responses,  while  the  dorsal 
hippocampus  is  more  involved  in  fear  learning  [55].  Because  the  sorting  of  behavioral  re¬ 
sponses  to  these  categories  is  not  feasible,  we  examined  the  whole  hippocampus  in  this  study. 
Since  the  goal  of  this  study  was  to  identify  QTGs  for  acute  stress  responsiveness,  the  assump¬ 
tion  was  that  steady-state  expression  of  DEGs  between  the  F344  and  WKY  strains  will  differ  in 
all  three  brain  regions  due  to  sequence  variations  in  the  candidate  QTGs. 

As  a  proof-of-concept,  we  hypothesized  that  the  expression  of  candidate  QTGs  in  the  brain 
of  the  congenic  strain  F344.WKY-SfresplO  would  parallel  that  of  the  WKY  and  differ  from  the 
F344  parental  strain.  This  pattern  would  mirror  the  strain  differences  in  behavioral  phenotypes 
obtained  either  in  the  DB  or  the  OFT.  We  identified  sequence  variations  within  12  genes,  of 
which  a  total  of  five  candidate  QTGs  were  identified.  Expression  of  three  of  the  genes  ( Cdkl4 , 
Drcl,  and  Gpatchll)  were  parallel  to  the  behavior  in  two  or  more  brain  regions,  while  the 
other  two  ( Cyplbl  and  Rbks )  were  parallel  to  the  behavior  in  at  least  one  brain  region.  Cdkl4 
belongs  to  the  cyclin-dependent  protein  kinase  family  and  is  responsible  for  cell  cycle  progres¬ 
sion,  including  transferase  activity,  transferring  phosphorous-containing  groups  and  protein 
tyrosine  kinase  activity  [56].  Drcl  encodes  a  key  component  of  the  nexin-dynein  complex  that 
regulates  the  assembly  of  ciliary  dynein  [57].  Gpatchll  is  involved  in  nucleic  acid  binding  [58]. 
Cyplbl  encodes  a  member  of  the  cytochrome  P450  superfamily  of  enzymes,  which  catalyze 
reactions  involved  in  drug  metabolism  and  synthesis  of  cholesterol,  steroids,  and  other  lipids 
[59-61].  In  retinal  endothelial  cells,  expression  of  Cyplbl  has  been  shown  to  reduce  intracellu¬ 
lar  oxidative  stress;  although  this  has  not  been  shown  in  neurons  [62].  Rbks  encodes  a  member 
of  the  carbohydrate  kinase  PkfB  family  and  is  known  to  catalyze  the  phosphorylation  of  ribose 
[63].  Although  none  of  these  genes  have  been  directly  implicated  for  acute  stress  responsive¬ 
ness,  their  genomic  location  within  QTLs  for  acute  stress  response  suggests  that  they  may 
underlie  some  common  mechanisms  of  these  phenotypic  traits. 

Interestingly,  the  individual  candidate  QTGs  are  link  to  stress-related  immunoregulatory 
genes,  including  115  (interleukin  5),  Btnl2  (butyrophilin-like  2),  Ifna2  (interferon  alpha  2),  and 
Ifnll  (interferon  lambda  1)  [64-67].  Both  115  and  Btnl2  can  activate  candidate  QTG  Drcl, 
which  is  known  to  encode  a  key  component  of  the  nexin-dynein  complex  that  regulates  the 
assembly  of  ciliary  dynein  [64,  65,  68].  While  Drcl  has  never  been  implicated  in  stress 
responses  or  stress-related  disorders,  115  has  been  reported  to  be  differentially  expressed  in  the 
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frontal  cortex  of  rats  exposed  to  acute  stress  [69].  In  a  human  study,  elevated  levels  of  7/5  were 
associated  with  an  increased  likelihood  of  major  depressive  disorder  [70, 71].  Additionally,  can¬ 
didate  QTG  Drcl  can  be  activated  by  Cd38  (Cd38  molecule),  which  encodes  a  multifunctional 
protein  involved  in  glucose-mediated  insulin  secretion  and  immune  system  functioning  [64,  65, 
72].  In  the  brain,  Cd38  regulates  the  secretion  of  the  neuropeptide  oxytocin  and  is  associated 
with  several  stress-related  phenotypes,  including  social  impairments  in  humans  such  as  autism 
spectrum  disorder  [73].  Additionally,  the  rs3796863  SNP  is  associated  with  heightened  stress 
sensitivity  and  predicting  social  anxiety  and  depression  in  humans  [73].  The  other  candidate 
QTG  Gpatchl  1  can  also  be  activated  by  Ifna2  and  Ifnll ,  which  are  both  interferon  immunosup- 
pressor  genes.  Interestingly,  Ifna2  is  a  pleiotropic  cytokine  that  triggers  immune  responses, 
hypothalamic-pituitary-adrenal  axis  abnormalities  and  disturbances  in  brain  metabolism  resem¬ 
bling  those  in  depressive  states  [74].  Ifna2  is  also  known  to  induce  memory,  concentration,  and 
mood  disturbances  when  administered  as  a  therapeutic  [74],  In  differentiating  neurons,  the 
expression  of  Ifna2  affects  their  response  to  inflammatory  cytokines,  which  is  consistent  with 
molecular  mechanisms  involved  in  schizophrenia  and  autism  spectrum  disorder  [75]. 

One  inherent  limitation  of  this  study  stems  from  the  heterozygous  region  of  the  congenic 
strain,  which  we  excluded  from  our  analyses.  Other  candidate  QTGs  could  be  mapped  to  this 
region.  Among  the  other  limitations  of  this  study  is  that  it  focuses  on  candidate  QTGs  with 
non-synonymous  sequence  variations  in  their  coding  regions.  It  is  known  that  the  large  major¬ 
ity  of  sequence  variations  are  in  non-coding  regions  that  may  act  as  cis-regulatory  and/or 
trans-acting  modules.  This  makes  identification  of  candidate  QTGs  for  behavioral  and  psychi¬ 
atric  phenotypes  more  difficult  [76].  However,  our  presumption  that  the  SNPs  within  the 
QTGs  affect  the  expression  of  these  genes  prompted  us  to  first  investigate  coding  region  non- 
synonymous  sequence  variations.  Furthermore,  among  inherited  gene  variations  in  humans, 
the  non-synonymous  SNPs  in  coding  regions  that  lead  to  changes  in  amino  acid  in  protein 
expression  are  most  relevant  to  human  inherited  diseases  [32],  Additionally,  these  candidate 
QTGs  with  cis-regulated  expression  changes  can  affect  gene  expression  in  trans,  as  described 
in  Bryois  et  al.  [77].  Future  studies  will  examine  sequence  variations  between  F344  and  WKY 
within  the  QTL  in  non-coding  sequences  conserved  across  species,  as  described  in  Yoshihara 
et  al.  [78].  We  will  also  focus  on  sequence  variations  in  microRNAs  with  known  gene  targets 
(www.targetscan.org  and  www.exiqon.com/microrna-target-prediction). 

Taken  together,  our  findings  indicate  that  strain  differences  in  acute  stress  responsiveness 
are  generalizable  across  multiple  behavioral  paradigms.  The  unique  approach  of  transcriptomics 
combined  with  sequence  variations  within  a  specific  QTL  in  the  different  parental  and  congenic 
strains  identified  unique  QTGs  that  might  contribute  to  variations  in  the  behavioral  responses 
to  acute  stress.  The  role  of  these  candidate  genes  in  the  behavioral  response  to  stress  should  be 
confirmed  in  future  studies  requiring  brain  region-  and  strain-specific  transgenic  approaches. 
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